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Abstract 

The geometrical coupling between strain and twist in double helices is investigated. Overwinding, 
where strain leads to further winding, is shown to be a universal property for helices, which are 
stretched along their longitudinal axis when the initial pitch angle is below the zero-twist angle 
(39.4°). Unwinding occurs at larger pitch angles. The zero-twist angle is the unique pitch angle at 
the point between overwinding and unwinding, and it is independent of the mechanical properties 
of the double helix. This suggests the existence of zero-twist structures, i.e. structures that display 
neither overwinding, nor unwinding under strain. Estimates of the overwinding of DNA, chromatin, 
and RNA are given. 



1 Introduction 

If one pulls a double helix structure by the end, one might think that it would unwind by the applied 
tension. In this paper we show why this is not always the case: A helix can unwind, overwind, or it 
can stay at its current twist (which we denote a zero-twist (ZT) structure). Overwinding is contrary 
to unwinding; unwinding is the de-twisting of the helices obtained by stretching the material. For the 
zero- twist structure there is no coupling from strain to twist. The existence of a twist-stretch coupling 
is a well-known phenomenon for helical steel wires [T] where it leads to unwinding, and design efforts go 
into designing rotation resistant wire rope when desired [2J [3] . 

The geometrical investigation presented below is based on the study of packed double helices modeled 
as two flexible tubes with hard walls. To be packed is defined by the constraint of the two tubes being 
in contact. Does this mean that the helices are stretched? No, generally not, stretching is one way to 
secure that a packed helix is obtained, however, for helices on the molecular size favorable molecular 
interactions can also make it more preferable to be packed than not. A detailed analysis of packed helices 
and their volume fractions showed that the helices with the highest volume fractions are noticeably 
similar to the molecular structure of DNA [4] ; this suggests that close-packing is at work as a structure 
forming principle. For the description of compact strings and tube models, the importance of one kind of 
optimum shape has been discussed by Gonzalez and Maddocks [5] and Maritan et al. [5] , and one related 
suggestion for the best packing of proteins and DNA has been considered by Stasiak and Maddocks 7 . 
A detailed analysis of the geometry of n-plies, and of their self-contacts, has been given by Neukirch and 
van der Heijden [8]. 

2 Model 

The close-packed (CP) structure with an optimized volume fraction has a pitch angle of 32.5° ]4i|: This 
structure that has a central channel is shown in Figure la. Under a pull, the pitch angle is increased 
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and the diameter of the central channel gets smaller, and eventually, the inner channel disappears at a 
pitch angle of 45° . Whether a helix overwinds or unwinds is then determined from the balance between 
the gain in length from the reduction in the helical radius versus untwisting. The crossing point - which 
we denote as the zero-twist angle - is at 39.4° (Figure lb) and is smaller than the 45°, where the helical 
radius becomes equal to the diameter of the tubes, and is maintained for all pitch angles above 45°. The 
45° motif, here denoted the tightly packed (TP) double helix, is shown in Figure lc. 



a) b) c) 




Figure 1: Different geometries of a double helix of tubes of fixed diameter D. a) Close-packed (CP) 
structure of pitch angle 32.5° measured from horizontal, b) Zero-twist (ZT) structure with a pitch angle 
of 39.4° . It is at the point between overwinding /unwinding, c) Tightly-packed (TP) structure of pitch 
angle 45°. Overwinding from stretching takes place from the a) to the b) confirmations; unwinding from 
b) to c). 

Geometrically, the double helix is given by two tubes of diameter D, whose centerline defines two 
helices with simple parametric equations. A helix is a curve of constant curvature, K, and torsion, t, and 
it can be specified by two parameters, for example a and H, where a is the helix radius (the radius of 
the cylinder hosting the helical lines) and H the helical pitch (the raise of the helix for each 2ir rotation). 
The tangent to each of the helical curves is at an angle v± (the pitch angle) with the horizontal axis, 
and it is determined by t&nv± = h/a, where h — H/2n is the reduced pitch. We say that the double 
helix is packed when the shortest distance between the centerline of one helical tube to the next one 
equals the diameter D of the tubes, i.e. the double helix is packed when the tubes are in contact. The 
volume fraction can be calculated using, as a reference volume, an enclosing cylinder of height H = 2nh 
and volume Ve = 2ir 2 h(a + D/2) 2 , and comparing it to the combined volume occupied by the two 
circumscribed helical tubes, Vh = ir 2 hD 2 / sinux- The volume fraction is the ratio of the two volumes, 
i.e. fv = Vh/Ve, which reads 

fy = 2(1 + ■ + I)" 2 (2.1) 

With this choice of reference volume the packing fraction depends only on the shape of the double helix 
structure, which can be described by one parameter, e.g. the pitch angle, v±. The maximum of fv defines 
the close-packed (CP) helix. For the double helix this maximum is at Vj_ — 32.5°, where f v = 0.796 
[4]. For the CP structure, the channel radius is about 17 % of a [4]. Generally, the radius of the central 
channel, which is given by Ri = a — D/2, is a decreasing function of v±; this can be seen from Figure 2 
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which shows 2a /D depending on the pitch angle. For v± > 45° there is no central channel as 2a /D = 1, 
see Figure 2. 




Figure 2: Graph showing the ratio 2a/ D as a function of pitch angle, v± [deg.], where a is the helix radius 
and D the diameter of the helical tubes. The tightly packed double helix has a pitch angle of vtp = 45°; 
it is the helix with the smallest pitch angle obeying the criterion that 2a = D. 



3 Results 

Consider a long straight segment of a double helix consisting of two long molecular strands each of length 
Lm- The length of the double helix is Hm = Lm smv± and the total twist is ®m — Lm cosv±/a. In 
Figure 3 the dimensionless ratio D9m /2Lm is shown as a function of the pitch angle. One can see that 
for vj_ < vzt there is overwinding while for vj_ > vzt there will be unwinding. We find numerically that 
v ZT = 39.4°. 

We can determine the amount of overwinding and unwinding in the following way. If a long double 
helical segment is stretched a bit, the pitch angle, v±, will change by a small amount dv±, and hence Hm 
changes by 

dHu — Lm cosuj_£foj_ (3-2) 

and Om by 

jr . T sin v. L M da 

a<OM = —Lm dv± ^-coswj^- — dv± (3-oJ 

a gt dv± 

so that 

d@M 1 1 da 

—— = — tan v ± - - — 3.4 
dtiM a a z av± 

If this derivative is positive, then the helix will overwind, and if it is negative, it will unwind. The 



derivative in Eq. (3.4) has dimension of inverse length. From a geometrical viewpoint it is more natural 
to look at the dimensionless function of v±, obtained by multiplying with the common radius of the 
tubes, (D/2), namely: 

DdQ ^ = - D t^ + ^( D ) (3.5) 
2 dH M 2a ± dv± \2a y ' 
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Figure 3: The total twist, 9m, for a long segment of the double helix; the dimensionless quantity D6m /2Lm 
is shown as a function of the pitch angle, vj_ [deg.J. The maximum value is obtained for the pitch angle 
vzt = 39.4° and mark the transition from overwinding to unwinding. At the ZT structure there is zero 
coupling between twist and strain. 



This equation can be given a simple interpretation. The first term is negative and determines the amount 
of unwind, while the second term is positive and determines the amount of overwind. The graph of this 
derivative, that dictates the coupling between strain and twist, is depicted in Figure 4. Notice that the CP 
double helix will always overwind since dQju /dHja > 0. At the close-packed structure, the overwind is 
(D /2)d®M /AHm = 0.665. The extension is therefore universally determined just by giving the diameter, 
D, of the tubes making up any close-packed double helix. At the zero-twist structure, vzt — 39.4°, 
there is neither overwinding, nor unwinding. For larger pitch angles the overwind, {D /2)dQM / oIHm, 
is negative and the double helix will unwind under strain. It is therefore crucial, that the pitch angle 
is below that of the zero-twist (39.4°) for overwinding to be observed, but it also indicates that elastic 
properties of the material are not essential to understanding the phenomenon. 



4 Discussion 



In the following we discuss some molecular examples. The phenomenon of overwinding in DNA was 
first observed in 2006, see Lionnet et al. [5] and Gore et al. [10] using magnetic tweezers to control the 
wringing [9] and optical tweezers to control the pulling [10] : For small deformations, DNA overwinds 
when stretched, i.e. it rotates counter to unwinding. During overwinding the extension of a long chain of 
DNA-B has been reported to be 0.42 ± 0.2 nm per 2ir rotation [5] and 0.5 nm per 2n rotation [TO]. Very 
recently, it has been suggested that in the absence of tension DNA is an order of magnitude softer [TT] , 

Using the above mathematical solution for the double helical structure of DNA we find the change of 
length AH to be determined by 

AH = S^A6 (4.6) 
d&M 

The diameter of the molecular tubes that make up the DNA helix is D = 1.15 nm, which is given from 
our previous analysis of the close-packed structures [3]. We then estimate AH per full 2ir turn to be 
7r(0.665) _1 x 1.15 nm = 5.4 nm, see Figure 4. Our result seems to support the findings of ref. [TT] , 
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Figure 4: Graph showing the calculated overwind of double helices (solid line), i.e. Eq. (3.5) as a function 
of vj_ [deg.J. A positive overwind means that the double helix will exhibit overwinding, while a negative 
overwind means that the double helix will exhibit unwinding. The zero-twist structure (ZT) is indicated 
with an arrow at vzt = 39. 4°, the close-packed structure (CP) is indicated by an arrow at vcp = 32.5°. 
The first derivative is discontinuous at vtp = 45° where the helix radius can not get smaller. The dashed 
line is the overwind for a triple helix, which has a zero-twist angle of 42.8° . 



The geometrical restriction imposed by base pairing and its influence on g?Om/^-Hm has not been taken 
into account. The numerical analysis has been performed for the symmetrical double helix where the 
close-packed structure has a pitch angle of 32.5°. The asymmetrical DNA-B has a close-packed pitch 
angle of 38.3° and, as one can show, a zero-twist angle of 41.8°. Theoretical work on understanding 
the overwinding of DNA has focused on constructing elastic models which show a negative twist-stretch 
coupling [12] and on incorporating stochastic effects [13] . One elastic model was considered by Gore et al. 
[TO] , and consists of a rod with a stiff helical wire (analogous to the sugar-phosphate backbone) attached 
to its surface. As this system is stretched, the inner rod decreases in diameter and the helix will overwind. 
Smith and Healey has argued that a linear material law is inadequate for the description and suggest a 
non-linear elastic rod |14) . 

For chromatin, the above results can be related to recent experiments in twisting chromatin fibers, 
see e.g. [T5] [IB]- For a close-packed 30 nm chromatin fiber, in the so-called two-start geometry, we 
estimate a tube diameter of 30/(2a/D + 1) nm= 30/(1.2 + 1) = 13.6 nm, where 2a/D is determined 
from Figure 2. For the close-packed 30 nm chromatin structure we then estimate AH per full 2tt turn 
to be 7r(0.665) -1 x 13.6 nm = 64 nm. It is interesting to note that the numbers reported in ref. [16] are 
measurements of AH for Xenopus chromatin per turn at a pulling force of 0.3 pN. Using the depicted 
data in ref. [16) we have estimated an average extension of ~ 60 ± 40 nm per turn. Here, we have 
assumed the two-start helix to behave like a tubular packed double helix - that is a view which ignores 
the intricate details of the structure, details which are discussed for example by Barbi et al. [25) . where 
elaborate mechanical models are described, including one which maintain its twist while being stretched. 

We have presented a simple geometrical explanation for overwinding of helices - an effect which has 
been observed before for the double helix of DNA and for chromatin, and which is contrary to usual 
unwinding. Our model of unwinding and overwinding can be applied to any symmetric double helix 
which is packed in the sense that the two helices touch each other, i.e. remain at the distance D from 
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each other. Packed double helical structures will show an overwinding behavior similar to those already 
observed, as long as their initial pitch angle is sufficiently small. Perhaps, the analysis will be relevant for 
other helical structures such as nanofabricated quartz cylinders [17] , fabricated twisted polymer nanofibers 
[T8] . and for the beautiful double helical structures formed from helical carbon nanotubes |19| . Further, 
the phenomenon may be important for some aspects of the working of molecular motors during gene 
expression and regulation |20j . The analysis presented in this paper is straightforwardly applicable to 
RNA double helices [21], which we therefore predict will show overwinding. Using a value of 26 A [22] for 
the molecular diameter of the double helix, we estimate an overwinding of 5.6 nm. Necturus chromatin 
fibers . 23] are known to pack as a double helix with a pitch angle of v± = 32 ±3° a value suggestive of being 
close-packed. Thus it follows that these chromatin double helices will overwind as well (other chromatin 
fibers with a different linker length would not necessarily overwind). Such predictions for overwinding 
and unwinding can nowadays be studied on single biomolecules using magnetic traps [22]. Furthermore, 
the derived geometrical expressions for overwinding are straightforwardly extended to helices with more 
than two strands. In Figure 4 we have shown the solution for a triple helix (dashed line) which has a 
zero-twist angle of 42.8°. 

Maybe one will even find examples, where Nature has build zero-twist structures, i.e. structures that 
display neither overwinding, nor unwinding. Chromatin with an appropriate linker length, and collagen 
are possible candidates for structures with such properties. 
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